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Abstract 



We consider pure glue QCD at (3 = 5.7, (3 = 6.0 and = 6.3. We 
evaluate the gluon propagator both in time at zero 3-momentum and 
in momentum space. From the former quantity we obtain evidence for 
a dynamically generated effective mass, which at (3 = 6.0 and (3 = 6.3 
increases with the time separation of the sources, in agreement with 
earlier results. The momentum space propagator G(k) provides fur- 
ther evidence for mass generation. In particular, at (3 = 6.0, for k < 1 
GeV, the propagator G(k) can be fit to a continuum formula proposed 
by Gribov and others, which contains a mass scale b, presumably re- 
lated to the hadronization mass scale. For higher momenta Gribov's 
model no longer provides a good fit, as G(k) tends rather to follow 
an inverse power law ~ l/k 2+1 . The results at (3 = 6.3 are consistent 
with those at (3 = 6.0, but only the high momentum region is acces- 
sible on this lattice. We find b in the range of three to four hundred 
MeV and 7 about 0.7. Fits to particle + ghost expressions are also 
possible, often resulting in low values for Xdofi Du t the parameters are 
very poorly determined. On the other hand, at (3 = 5.7 (where we can 
only study momenta up to 1 GeV) G(k) is best fit to a simple massive 
boson propagator with mass m. We argue that such a discrepancy 
may be related to a lack of scaling for low momenta at (3 = 5.7. 

From our results, the study of correlation functions in momentum 
space looks promising, especially because the data points in Fourier 
space turn out to be much less correlated than in real space. 



1 Introduction 

The possibility of studying nonperturbatively on the lattice gauge-dependent 
quantities provides in principle a unique tool to test QCD at the level of the 
basic fields entering the continuum Lagrangian. From this point of view, the 
gluon propagator in the pure glue theory is perhaps the simplest quantity. 
From its study one expects to obtain among other things a better under- 
standing of the infrared behavior of the theory and of the mechanism of 
gluon confinement. As a result we may also hope to acquire a better un- 
derstanding of the hadronization phenomena and of the glueball spectrum 
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Let us first review what is known from perturbation theory: consider the 
expression (in Minkowski space-time) 

Df v {k) = -iJ d'x <0| T[A;(x) 4(0)] |0) e*- (1) 

From the Faddeev- Popov quantization in a class of covariant gauges one gets 
the simple Slavnov identity 

Wk v Df u {k) = -ia 5 ab (2) 

where a is the gauge parameter. From Lorentz covariance, the general solu- 
tion to (|2|) can be written as 



D%{k) = -i5 ab 



k^k v \ 1 k^ky 

g»v - -j-r j 1 + n ( fc2ja ) + a ~gr 



F (3) 



Eq. (0) implies that the longitudinal part of the propagator gets trivially 
renormalized, so that the vacuum polarization IT(/c 2 , a) just renormalizes the 
transverse propagator. The renormalization constant Z 3 for the gauge fields, 
Ai = Zl /2 Af is defined as 

Z 8 - 1 = l + II(£,a) (4) 

where A is an ultraviolet cut-off and k 2 = — /i 2 is a spacelike value of the 
momentum. In general Z% is gauge dependent. 

One can then rewrite the unrenormalized transverse gluon propagator as 



Dl ab {k 2 )\ k , = _^ = ± Zi L, + ^ ) 6°» (5) 



Unlike QED, in QCD we know from the properties of the /3-function that 
perturbation theory can only be applied to the study of the large-fc 2 behavior 
of Green's functions. In such a region one gets at the 1-loop level M 



l + jLfll-a). 3W-) (6) 
16vr 2 V3 J V /i 



In general, at each order in g 2 perturbative corrections depend logarithmi- 
cally in k, so that in the deep Euclidean region (k 2 — > oo) gluon propagators 
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behave essentially in the same way as the photon propagators, as long as one 
considers a finite order in perturbation theory. 

In the infrared region things are probably different. In fact, in the simpler 
QED case perturbation theory is still reliable, and one can evaluate order 
by order the vacuum polarization function. In particular, one finds that 
UP ED (k 2 = 0) is finite if all the fermions are massive (we omit from now on 
the explicit a— dependence of vacuum polarization functions). This implies 
that the k 2 = pole of the free photon propagator is still present after 
radiative corrections are taken into account, so that the photon remains 
a massless particle. Such corrections only affect the residue at the pole, 
resulting in charge renormalization. 

On the other hand, consider QCD without quarks. The corresponding 
Lagrangian does not contain any mass scale in 4 dimensions, yet if the theory 
is confining a mass scale must be dynamically generated in some way, since 
the confinement potential V(r) = Kr contains such a scale. Such a mass M 
cannot be generated in perturbation theory, since it must satisfy 



M(g, fi) = fi exp 



a dg 



where /x is a renormalization scale /x. For small g, one has 

const 



(7) 



M(jg,fi) ~ n exp 



9 2 



when g — > (8) 



so that M(g) has an essential singularity at g = 0. 

Such a mass scale may show up in the vacuum polarization function for 
the gluon. Indeed, while contributions to U(k 2 ) proportional to finite powers 
of g have for large k 2 a logarithmic momentum dependence, non-perturbative 
effects may generate terms in U(k 2 ) proportional to negative powers of the 
squared 4-momentum, for instance 

n^.-^+^+o^ ( 9) 

where m(g,fi), b(g,p) have the dimension of a mass and depend nonanalyt- 
ically on g. For example, the case 6 = 0, m/fl gives rise to a mass pole in 
the gluon propagator and corresponds to the standard Schwinger mechanism 
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In the above formula 0(g 2 ) denotes contributions which can be repre- 
sented as a power series in g 2 . Such series, when truncated to a finite order 
in g 2 , behave like polynomials in ln(k 2 /n 2 ) for large k 2 . On the other hand, 
the sum of the contributions to all orders in g 2 can generate an anomalous di- 
mension 7 [f|] ; in such for k 2 — > oo the propagator behaves like l/k 2+1 . 
The nonperturbative behavior of the Euclidean gluon propagator has been 
investigated in the continuum by many authors with different methods and 
in different gauges 0, ||, [7|, [8|, |9], |10| . In some of these attempts a very singular 
gluon propagator was found, behaving like k~ 4 in the limit k 2 — * fL0| , and 



a confining property was inferred from such behavior. On the other hand, 
other work points towards the elimination at the non-perturbative level of 
the singularity at k 2 = of the propagator || [7|, [|, as a consequence of 
dynamical mass generation. 

In particular, a very peculiar momentum dependence, consistent with 
the above scenario, has been predicted as arising from a modification of 
the standard path integral Faddeev-Popov formula in the Landau gauge by 
the introduction of a nonperturbative gauge- fixing procedure p, |[. Such 
improved implementation of the Landau gauge is expressed by the equations 

d ■ A = and FP[A) > (10) 

where FP[A] is the Faddeev-Popov operator in the Landau gauge, which in 
general is not positive definite. The positivity requirement in (|T0"| ) can be 



seen as a recipe to get rid (although not completely |TTJ) of Gribov copies 
[|]. In the gauge (|T0|) , the (transverse) gluon propagator in momentum space 
has been argued to be of the form || 

where b is a dynamically generated mass scale. Eq. (|TT|) corresponds to the 
case m = 0, b ^ in (Q). It implies that in the continuum 

G(k = 0,t) « e"75* \cos{^=t) - sin(^=t)^j (12) 

Remarkably, the same predictions were also obtained in the study of Schwinger- 
Dyson equations 0. 
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The above expression lends itself to intriguing speculations: the absence 
of any particle singularity on the real k 2 axis predicts the absence of an 
asymptotic gluon state. It may describe a short-lived excitation, giving rise 
to a gluon jet. In this framework, the mass scale b appearing in the above 
formulas may perhaps be interpreted as a hadronization scale. 

The lattice gluon field can be defined as |jl2|: 



Ar(n) ^M-Uli.n)_, T IUM-U>)\ 



2ia 3 \ 2ia 



where a is the lattice spacing. Thus the lattice gluon propagator in x-space 
is the expectation value of: 

G^(x,y)=Tr{A^x)A u (y)) (14) 

An important point is that on the lattice one can define and implement 
the analogue of the gauge condition (|10|). In fact, given any link configuration 
{£/}, one can define a function of the gauge transformations g on {U} 

Fu[g\ = ~T. Re Tr (u°(n) + Uf{n - £)) , (15) 

where V is the lattice volume and U 9 indicates the gauge-transformed link 
U 9 (n) = g(n)U^{n)g\n + ft). An iterative minimization of Fjj[g\ obtained 
by performing suitable gauge transformations generates a configuration {U 9 } 
which satisfies the lattice version of (|10D, defined in terms of a lattice Faddeev- 
Popov operator 0|. This is just the Hessian matrix associated with Fu[g], 



that is, j Fu J 9 ^ . As we have already mentioned for the continuum, such a 
gauge is not completely free of Gribov copies. This is also true on the lattice, 
and corresponds to the fact that for a fixed configuration {U} the function 
( |15D may have several local minima (see, for example, [p~f ). 

In general, Zwanziger [^] showed that on the lattice one has qualitatively 
the same scenario for the Landau gauge as in the continuum. Indeed, there 
exists a bounded region Q, defined by the positivity requirement for the lat- 
tice Faddeev-Popov operator, which satisfies bounds analogous to the ones 
derived for the continuum model. Considering then a restriction of the func- 
tional integration to the region Q, Zwanziger was able to obtain predictions 
for the lattice gluon propagator consistent with the continuum ones given in 
(0) and (0). 
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At this point it is natural to try and test numerically predictions like (O) 
and (H). 

Numerical studies have been performed in the past years for the zero spa- 
tial momentum Fourier transform of (|14D, namely G{k — 0, t) = J2i=i Gaik 



I). I ) |TJ, [15|, [l(J. These studies reported evidence of an effective gluon mass 
that increases with the time separation for short time intervals. This feature, 
which would be unacceptable for the propagator of a real physical particle 
since it violates the Kallen-Lehmann representation, is in qualitative agree- 
ment with the continuum prediction (|i~2l) and may be in principle acceptable 
for a confined particle 0, |T2fl . Another lattice approach to the gluon mass 
was given in [17] . 



Our work aims to test at a more quantitative level continuum predictions 
and to extend the above results through the study of the gluon propagator 
at nonzero momenta. 



2 Numerical Results 

We study pure glue QCD on 16 3 x 40 and 24 3 x 40 lattices at (3 = 6.0, on a 
24 4 lattice at /3 = 6.3 and on a 16 3 x 24 lattice at f3 = 5.7. 

2.1 Technical Remarks 

It is worth remarking that, unlike simulations involving quenched quark prop- 
agators, evaluations of purely gluonic correlation functions can take full ad- 
vantage of the translational symmetry of the theory in order to improve 
statistics. On the other hand, such quantities turn out to be very sensi- 
tive to the numerical accuracy of gauge fixing. Empirically, at (3 = 6.0 and 
(3 = 6.3 we find that when the minimization of Fu[g] has reached an ac- 
curacy such that in O(50) iterations Fjj[g] changes less than « .05%, then 
the signal for the propagators is sufficiently stable against additional gauge 
fixing. In other words, the variation in each data point for the propagator 
arising from additional gauge fixing is typically much smaller than the final 
error bar associated to the data point. On the contrary, we will see that at 
P = 5.7 our stability requirement for F v [g\ does not suffice to guarantee a 
completely stable propagator. 

Even at (3 = 6.0 and (3 = 6.3, though, when our empirical criterion gets 
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satisfied the system has not yet reached complete equilibrium, in spite of the 
fact that such an accuracy is roughly one order of magnitude better than the 
typical one adopted in simulations of hadron phenomenology. This can be 



seen by performing the following test [12]: in the Landau gauge On An = it 



follows from the periodic boundary conditions that 

A (t) = J2 A)0M) (16) 

X 

should not depend on t. In lattice language, this means that once the Landau 
gauge has been numerically implemented in a configuration, then the sum 
over the sites in a fixed timeslice of the time component of the gauge field 
should be the same on each timeslice. In Fig. [I] we plot one of the diagonal 
elements of the matrix A (t) as a function of t for one of our gauge-fixed 
configurations on the 16 3 x 40 lattice at /3 — 6.0. This test shows that 
even when the accuracy of our numerical gauge fixing is sufficient for the 
gluon propagator, there are other quantities for which the gauge fixing need 
not be adequate. As a consequence, it would be dangerous in a calculation 
to rely on some standard a priori criterion when estimating the required 
precision for the gauge fixing, since such precision strongly depends on the 
specific observable under consideration. We remark that the gauge-fixing 
test provided by A (t) is in fact more stringent than the one used in \T2\. 



2.2 Results at (3 = 6.0 and /3 = 6.3 

Here we give results for a set of 25 configurations of size 16 3 x 40 at (3 = 6.0, 
a set of 8 configurations of size 24 3 x 40 at (3 = 6.0, and finally a set of 20 
configurations of size 24 4 at (3 — 6.3. 

As a first step we have evaluated G(k = 0, £); our results confirm that 
this propagator exhibits a massive decay in time, with an effective mass 
am(t) = ln(G(k = 0, t)/G(k = 0, t + a)) that increases with t for short times. 
In Fig. we plot am(t) versus t with jackknife errors [18] on the 16 3 x 40 
lattice. Assuming the value of the inverse lattice spacing or 1 = 2.1 GeV at 



(3 = 6.0 Ul9[ , the effective gluon mass m(t) ranges approximately between 
220 and 870 MeV. A similar behavior is observed for the 24 3 x 40 lattice, 
with the effective mass ranging between 240 and 1200 MeV. 

At this point we have attempted to fit G(k = 0, t) to the continuum form 
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( |i2|) and to the form commonly referred to as particle + ghost, that is 

G(k = 0, t) = G x exp(-Mit) + C 2 exp(-M 2 t) (17) 

where C 2 is constrained to be negative. 

As is well known, the data points obtained from a Monte Carlo simulation 
are in general statistically correlated; in the present case, the correlated data 
are the values of the propagator G(k = 0,t) at different timeslices. For this 
reason, one should perform y 2 fits by using the definition of x 2 which involves 



the full covariance matrix |20 



By inspection of the covariance matrix for G(k = 0,t), it turns out that 
the off-diagonal matrix elements are typically of the same size as the diagonal 
ones, which means that our data points are highly correlated in t. Conse- 
quently, x 2 fits are not well controlled because the covariance matrix is nearly 
singular. Much higher statistics would be required to get well-behaved fits. 

As it is not possible in this case to make fits using the full covariance 
matrix, we are forced to use the naive definition of x 2 \ where the data points 
are simply weighted by their standard error bar. Of course, in such an 
approach x 2 h as no simple relation to "goodness of fit." 

In this approximation, it turns out that Gribov's formula provides, for 
small t, a good fit to G(k = 0,t), better than the one obtained from the 
particle + ghost expression. In fact, although Xdof is not a reliable indicator 
of the goodness of a fit when the covariance matrix is not taken into account, 
nonetheless the relative Xdof °f two fits provides some indication of which fit 
is better. On the 16 3 x40 lattice at (3 = 6.0, the lowest value for x\of obtained 
from a fit to Gribov's formula is Xdof = • 18, while from the 4-parameter fit 
to particle + ghost one gets Xdof = -34- Moreover, the latter kind of fit is 
much less stable against varying the initial guess for the fit parameters, i.e. 
many local minima for x 2 can be found. Using a -1 = 2.1 GeV, we obtain 
for the b parameter in (|12D b = 237 ± 7 MeV, Xdof = -18; where the error on 
b is a jackknife one. Of course one does not get good agreement by using a 
conventional 4-parameter double exponential form, that is if one constrains 
C 2 in ( |T7D to take positive values. Indeed, in this case the effective mass 
would always decrease with t, in contrast to what is observed. 

The statistical difficulties mentioned above forbid a complete analysis of 
G(k = 0,t); in particular, we cannot effectively study the large-t region, 
where according to (|T2]) the propagator may become negative and then os- 
cillate. 
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Our analysis of the momentum space propagator G(k) = X^ =1 

does not have the same statistical difficulties as in the case of G(k = 0,t). 
This quantity is obtained by performing explicitly the lattice Fourier trans- 
form of the propagator (|Hp. 

In fact, the covariance matrix associated with G(k) turns out to be much 
more "diagonal" than the one for G(k = 0, £); in other words, the data points 
are much less correlated in momentum space than they are in t. We find that 
G(k) is very well determined in a significant interval of physical momenta, 
ranging from the lattice infrared cutoff in the time direction k = 2ir/N t a 
(N t = 40 at (3 = 6.0) up to k pa 2 GeV, which is in fact roughly the value of 
the ultraviolet cutoff a -1 (see Fig. As a consequence we have been able to 
obtain good fits to G(k) taking into account correlations. Fig. |3| also shows 
that the data points from the two lattices at (3 — 6.0 are in good agreement 
over most of the momentum range. 

At this point we attempt to fit G(k) to the continuum formula ( p"T|) and, 
for a comparison, to a standard massive propagator G mass {k) = A/{k 2 + m 2 ). 
The fits are performed including in \ 2 the covariance matrix. Consider first 
the 16 3 x 40 lattice: it turns out that for momenta up to ~ 1 GeV formula 
( |TTD fits the data well (see Fig. Indeed, for such a fit we obtain Xdof = 1-3 
and b = 341 ± 12 MeV, assuming a -1 = 2.1 GeV. The error is given in 
terms of the parameter covariance matrix [20|, and we are not including the 
uncertainty in the value of a -1 . 

We compare this result to the best fit that one can obtain from the stan- 
dard massive propagator, for which we obtain Xdof = 2-9 Q. 

On the other hand, in the range of momenta 1 — 2 GeV the formula (|Tl|) 
and the massive propagator expression fit the data very poorly, resulting in 
Xdof > 10. A good fit in this range is obtained by assuming an inverse power 
law behavior G(k) = A/k 2+ "> (see Fig. §). We obtain 7 = .7 ± .2, Xdof = L 
The fact that ([TT|) does not fit the high momentum region well presumably 
explains the discrepancy between the value of b obtained from G(k) at low 
momenta and that obtained from G(k = 0,t), which depends on the complete 
range of momenta. 

The situation is qualitatively similar on the 24 3 x 40 lattice. However, 



*We have also attempted a fit to a form proposed in ||. Such a form in the low 
momentum region is very similar to a standard massive propagator. The standard massive 
form appears to provide a better fit. 
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presumably due to the lower statistics, there is clear dispersion in the data 
(see Fig. |3|), which means that no fit will be particularly good. In fact, up to 
1 GeV the best fit is given by (JTT|) (see Fig. and we obtain b = 333 ± 12 
MeV, Xdof = 5-8. Between 1 and 2 GeV we are unable to perform fits with 
the covariance matrix, due to the poor statistics. Without the covariance 
matrix, an inverse power law with 7 « .6 again reproduces the data well (see 
Fig. |). 

Summarizing, the data for G(k) at (3 = 6.0 on the 16 3 x 40 lattice up to 
momenta of order 1 GeV prefer somewhat formula (|TT|) , which describes the 
mass generation a la Gribov, over a standard massive propagator. In the 
momentum range 1-2 GeV the propagator is best reproduced by the inverse 
power law behavior G(k) m l/k 2+1 , where 7 could be interpreted as the 
anomalous dimension of G(k). Indeed, our results in such a region could also 
be sensitive to the lattice ultraviolet cutoff a -1 . The fact that the (3 = 6.3 
results (described below) agree with those at (3 = 6.0 gives some indication 
that such cutoff effects are not overwhelmingly large. An additional test for 
lattice artifacts is described at the end of the next section. 

On the other hand, the results from the 24 3 x 40 lattice, due to the poor 
statistics, do not provide by themselves a strong confirmation for the above 
picture. Still, in the low momentum region the Gribov fit works better than 
a standard massive one and gives a value for b which is consistent with the 
one from the 16 3 x 40 lattice, but we get a high value for Xdof- Moreover, in 
the higher momentum region we are unable to use the covariance matrix. 

We consider now a set of 20 configurations on a 24 4 lattice at (3 = 6.3. 

Starting again from the evaluation of G(k = 0,t), we give in Fig. |6| the 
effective gluon mass with jackknife errors. Again, the effective mass am(t) 
appears to increase with t. Assuming a -1 = 3.2 GeV at (3 = 6.3 m{t) 
ranges approximately between 200 and 670 MeV. 

It turns out that we have the same statistics problem as at (3 = 6.0: the 
data for G(k = 0,t) are highly correlated in t, so that it is impractical to 
make fits using the full covariance matrix. 

Next we consider the propagator G(k) in momentum space. We again 
find that correlations between different values of k are much smaller than 
the correlations in t, so that the covariance matrix is now well behaved. The 
best fit to G(k) is obtained from the inverse power law G(k) = A/k 2+1 (see 
Fig. 0). We assume a~ l = 3.2 GeV U and get 7 = .68 ± .08, xLf = L9 - 

Given the momentum range of such a fit, which starts around 1 GeV and 
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goes up to momenta of the order a -1 , such a result is quite consistent with 
the behavior observed at f3 — 6.0. In particular, the value for 7 agrees with 
the one obtained on the 16 3 x 40 lattice. 

Both at (3 = 6.0 and f3 = 6.3 we have also attempted 4-parameter fits, in 
terms of the particle + ghost form G(k) = A/(k 2 + m 2 ) + B/(k 2 + m 2 ,), with 
A > and B < 0. Although low values of Xdof can sometimes be obtained, 
the determination of the fit parameters is very poor. 

2.3 Results at f3 = 5.7 and Comparison of Lattices 

Here we give results for a set of 16 configurations on a 16 3 x 24 lattice at 
(3 = 5.7. After implementing our standard level of gauge fixing, we find 
that the data for G(k = 0, t) are, as usual, not suitable for \ 2 fits with the 
full covariance matrix. When fit without the covariance matrix, the data 
seem to be best reproduced by a simple decreasing exponential, unlike what 
happened on the other data sets at weaker couplings. Correspondingly, there 
is no clear evidence for an increase of the effective gluon mass am(t) on a 
significant time interval (see Fig. |8|). 

In addition, the results for G(k) on this lattice differ from what we ob- 
served at weaker couplings. Assuming a -1 = 1.2 GeV, our data for G(k) cover 
a momentum range up to roughly 1 GeV. The best fit to such data (including 
the full covariance matrix) is provided by a free massive boson propagator 
G(k) « 1/k 2 + m 2 . We obtain m = 590 ± 30 MeV, with xLf = 1-4. 

Since from the above analysis the results at (3 — 5.7 seem to differ signif- 
icantly from those at (3 = 6.0 and (3 = 6.3, we have investigated the effect of 
improving the accuracy of the gauge fixing. After 300 additional sweeps on 
each configuration, we repeat our measurements. 

For G(k = 0, t) it turns out that while the central values of the individual 
data points do change significantly, the pattern of the effective mass am(t) 
is basically the same. Not surprisingly, the best fit is still provided by a 
simple exponential, with a value for the mass parameter which is close to 
the one obtained before the additional sweeps; in other words, to a first 
approximation the effect of additional gauge fixing has been a rescaling of 
the data points. As far as G(k) is concerned, the individual data points are 
not much affected by the additional gauge fixing. In fact, the best fit to the 
data is still provided by a free massive boson, but now we obtain m = 630 ±40 
MeV, with Xdof = 3-3. Moreover, the error bars of some data points have 
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become slightly bigger (see Fig. g). 

Summarizing, we find that the qualitative behavior of the data is not 
changed by the additional gauge fixing, maintaining the discrepancy with 
the other lattices. However, since the value of Xdof m the ^ for G{k) has 
increased as well as some individual error bars for the data points, the ad- 
ditional gauge fixing seems to have increased the "noise" in the evaluation 
of G{k). This may be related to the existence of Gribov copies in our gauge 

0- 

To go back to the main issue, which is the discrepancy between G(k) at 
P = 5.7 and at weaker couplings, our present data give some indication that 
such a discrepancy may be due to a lack of scaling at (3 = 5.7. To illustrate 
this point we have plotted for each data set the quantity 100 x G(k)/G(k = 0) . 
Such a quantity should not depend on /3, as long as one is in the scaling region. 
On the other hand, when comparing G(k)/G(k = 0) at different (3 one has 
an uncertainty related to the determination of the horizontal scale, i.e. the 
scale of physical momenta, depending on the values of a~ l . With this caveat 
in mind, we show in Figs. [10] and [TTJ the quantity 100 x G(k)/G(k = 0) as 
obtained from our different data sets. (For the lattice at (3 = 5.7 we use the 
data after the additional gauge fixing.) We again assume a -1 = 1.2 GeV at 
(3 = 5.7, a" 1 = 2.1 GeV at (3 = 6.0 and a' 1 = 3.2 GeV at (3 = 6.3. 

From the above mentioned figures it turns out that in general the data 
agree well, but significant deviations are observed in the low momentum 
region; in particular, the data at (3 = 5.7 for momenta up to m 650 MeV are 
rather different from the corresponding data at (3 = 6.0 (see Fig. [10]). This 
may account for the observed discrepancies, although further analysis and 
better statistics are needed to clarify this issue. 

At this point it is worth mentioning that we have checked the stability 
of our fits for G(k) versus the use of continuum and lattice formulas. In 
fact, especially when the the data points correspond to lattice momenta of 
order a -1 , one may wonder whether more accurate fits could be obtained by 
using lattice versions of our formulas. One can devise such expressions by 
substituting k 2 with 2/ a 2 J2n (1 — cos(k^a)) in the continuum ones. 

We have also checked the stability of the results under change from co- 
variant to noncovariant fits. It turns out that when the covariance matrix is 
well behaved, which is what typically happens on our lattices for G(k), the 
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fit is not very sensitive to the different definitions for x 2 - Moreover, the fits 
are not sensitive to the difference between continuum and lattice formulas, 
not even in the higher momentum region. This is shown in Fig. O for the 
24 4 lattice at (5 = 6.3, where we compare lattice and continuum versions of 
the best fit to an inverse power law, with and without the full covariance 
matrix. 

One gets from the fit to the lattice formula, without covariance matrix, 
7 ~ .69, and from the fit to the continuum formula, without the covariance 
matrix, 7 m .61. These numbers should be compared to 7 = .68 ± .08, which 
is the result we obtained from the fit to the continuum power law which 
included the full covariance matrix. 



3 Conclusions 

Let us summarize our results. 

We have evaluated G(k = 0,t) and G(k) on four different lattices at 
three different values of (3. All the data provide evidence for dynamical mass 
generation, in agreement with previous results for G(k = 0,t). 

1. For G(k = 0,t) the data at /3 — 6.0 and (3 = 6.3 show an effective 
gluon mass which increases with the time separation, for short times, 
in agreement with previous results. At (3 = 5.7 there is no clear sig- 
nal for an increasing effective mass over more than 1 time slice. A 
more detailed analysis of G{k = 0,t) is not possible, due to statistical 
difficulties. 

2. The most interesting quantity is the momentum space propagator G(k), 
which we have evaluated for the first time. It is characterized by a clear 
numerical signal; in particular, it turns out that the data in momentum 
space are much less correlated than in real space. At f3 = 6.0 the 
propagator G{k) up to momenta « 1 GeV can best be described by 
the dynamical mass generation mechanism a la Gribov; however, a 
description by a standard massive propagator is not ruled out. At 
higher momenta the propagator is best described by an inverse power 
law, which could be interpreted in terms of an anomalous dimension. 
The results at (3 = 6.3 are consistent, but only the high momentum 
region is accessible because of the small physical size of the lattice. 
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Thus the behavior of the propagator for the above values of (3 could be 
summarized by the determination of a mass scale b and an anomalous 
dimension 7. In spite of the uncertainties in the determination of a -1 , 
our numerical values for 7 are consistent between (3 = 6.0 and (3 = 6.3 
lattices, when assuming a -1 =2.1 GeV at /3 — 6.0 and a -1 = 3.2 
GeV at P — 6.3. For b, the results are consistent between the two 
different volumes at /3 — 6.0, which are the only lattices on which it is 
determined. 

On the other hand, at (3 — 5.7 G(k) is best fit, up to momenta 1 GeV, 
by a simple massive propagator. We argue that such a discrepancy may 
be related to a lack of scaling between (3 = 5.7 and 6.0 for G(k) at low 
momenta. 

In some cases, 4-parameter fits to particle + ghost expressions for G(k) 
result in small values of Xdofi but as the parameters are typically very 
poorly determined, such fits are not very illuminating. 

In conclusion, the lattice study of gluon correlation functions, in the 
gauge-fixed framework which we have discussed, appears to be technically 
possible, although challenging. In particular, the feasibility of the study in 
momentum space is very promising for future applications. 

A very careful analysis of systematic lattice effects and statistical errors 
is necessary. At the same time, a study of the gauge dependence of the mass 
scales related to the propagator is also in order. 
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Figure Captions 



Figure 1: A (t) vs. t for one gauge- fixed configuration on the 16 3 x 40 lattice 
at = 6.0. 



Figure 2: Effective gluon mass in lattice units vs. t on the 16 3 x 40 lattice 
at = 6.0. 



Figure 3: Momentum space propagator vs. k in GeV on the 16 3 x 40 and 
24 3 x 40 lattices at p = 6.0. We assume a" 1 = 2.1 GeV. 
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Figure 4: G(k) vs. k in GeV and fit to the form ( |lT|) on the 16 3 x 40 and 
24 3 x 40 lattices at (3 = 6.0. 



Figure 5: G(k) vs. k in GeV and fit to an inverse power law on the 16 3 x 40 
and 24 3 x 40 lattices at (3 = 6.0. In the latter case, the covariance matrix is 
not included in the fit. 



Figure 6: Effective gluon mass in lattice units vs. t at /3 = 6.3. 



Figure 7: The best fit of G(k) vs. k in GeV to an inverse power law (solid 
line) on the 24 4 lattice at (3 = 6.3 



Figure 8: Effective gluon mass in lattice units vs. t on the 16 3 x 24 lattice 
at (3 = 5.7. 



Figure 9: G(k) vs. k in GeV and fit to a massive boson propagator before 
and after additional gauge fixing on the 16 3 x 24 lattice at (3 — 5.7 



Figure 10: 100 x G(k)/G(k = 0) on the different lattices (low momentum 
region) . 



Figure 11: Same as Fig. [H], in a higher momentum region. 

Figure 12: G{k) vs. k in GeV on the 24 4 lattice at (3 — 6.3 and the best fits 
to an inverse power law 
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